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Abstract: In classical statistics, the population mean is estimated using determinate, crisp data value when 
auxiliary information is known. These estimates can often be biased. The main objective of this study is to 
introduce the neutrosophic estimator with the minimum mean squared error (MSE) for the unknown value of 
the population mean as well as overcome the limitations of classical statistics when dealing with ambiguous 
or indeterminate data. Neutrosophic statistics was introduced by Florentin Smarandache. It is a generalisation 
of classical statistics that addresses ambiguous, unclear, vague, and indeterminate data. In this study, we have 
proposed neutrosophic almost unbiased estimator that use known neutrosophic auxiliary parameters to 
estimate the neutrosophic population mean of the primary variable. Equations for bias and mean squared error 
are calculated for the suggested estimators up to the first order of approximation. The proposed estimator 
performs better than the other existing estimators with respect to the MSE and percent relative efficiency 
(PRE) criteria. The estimator with the highest PRE or lowest MSE is advised for practical utility in various 
kinds of application areas. The theoretical conclusions are validated by the empirical analysis, which made 
use of the real data sets. 


Keywords: Classical statistics, neutrosophic statistics, auxiliary information, bias, mean squared error, percent 
relative efficiency. 


1. Introduction 


One of the purpose of sampling theory is to estimate study variable on the basis of auxiliary 
information with more precision. In sampling theory, we estimate population parameter on the basis 
of auxiliary variable that has a correlation with study variable. Some of the popular methods used 
to obtain these objectives to estimate population parameter using ratio, product, regression method 
of estimation. 


Ratio and product estimators are frequently used in survey sampling for getting an improved 
estimator of population mean given that the auxiliary information is known. The concept of ratio 
estimator was introduced by [1] as a method for utilising auxiliary information. Ratio estimator 
compute ratios between the sample means of the study variable and the auxiliary variable and using 
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known auxiliary data on a variable that is linearly related to the variable under study, is used to 
estimate the population mean. 


Whereas, the product estimator, an additional technique for adding auxiliary information into 
estimation, was put up by [2]. In the product method of estimation, the sample means of auxiliary 
variable and study variable are multiplied, and the result is used as the foundation for estimation. 
Over time, several improvements have been made to the classical ratio and product estimators. 


One improvement was proposed by [3], they proposed exponential ratio and product type estimators 
for estimating unknown population mean of the study variable using auxiliary information. 


Several researchers such as [1,5,6,7] accomplished some outstanding work in this area. However, 
our primary focus here is on proposing almost unbiased estimator using neutrosophic information. 


Classical statistics deals with data form by crisp numbers and there is no uncertainty in 
measurement of the observation. But there are so many cases when indeterminacy, vagueness, 
ambiguousness exists in the data. In day-to-day life where we have indeterminate data, for 
illustration, the measurement of temperature anywhere, decision support systems, information 
systems, the weather forecast, economy growth and more. In such cases we need to introduce new 
techniques to handle with such type of indeterminate data. Fuzzy logic is one method for resolving 
these kinds of issues regarding ambiguous, imprecise or uncertain observation but still, it ignores 
indeterminacy. In this situation, neutrosophic techniques are more suitable which address both 
randomness and indeterminacy. 


While classical statistics deals with determinate and crisp values, neutrosophic statistics deals with 
sets values. As an extension of classical statistics, neutrosophic statistics is identical as classical 
statistics in the absence of indeterminacy, [9]. The evaluation of both determinate and indeterminate 
components of the observations is possible with neutrosophic logic, which is an extension of fuzzy 
logic. It is particularly useful for analysing data that is imprecise or unclear, as noted by [9]. 


Neutrosophic statistics can be used to solve a number of real-world issues as neutrosophic statistics 
to the analysis of road traffic accidents, [10]. [11] introduced a goodness-of-fit test within the 
framework of neutrosophic statistics and [12] concentrated on using neutrosophic analysis of 
variance on data from university students. In neutrosophic statistics, the sample size may not be 
precisely determined and not known as the exact number, [13]. Methods based on fuzzy logic are 
rapidly evolving and frequently used in decision-making settings, [14]. Further advancement of 
fuzzy sets is complex fuzzy sets, and a complex neutrosophic set is their generalized form. 
Estimation using neutrosophic techniques are a relatively new and mostly unexplored area 
compared to estimation challenges in classical probability sampling designs with definite and well- 
defined data, [15]. 


Classical statistics struggled to handle indeterminacies that arise from ambiguous, imprecise or 
incomplete frequency-related data, because it involved values that were uncertain and could include 
a range of possible measurements. Therefore, neutrosophic statistics is used as a more adaptable 
and extension of classical statistics in uncertain or indeterminacy situations. However, there has 
been limited attention given to neutrosophic sampling as compare to classical statistics for 
estimating the population parameters. 


Frame of our paper is constructed in the following manner: Section 1 mentioned Introduction. 
Section 2 discussed about observation and terminology of neutrosophic statistics. Section 3 outlines 
some existing neutrosophic estimators. Section 4 presents the proposed neutrosophic almost 
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unbiased estimator. Section 5 mentioned some member of the proposed estimator. Section 6 
presents empirical study. Section 7 is devoted to discussion whereas section 8 discusses the 
conclusion. 


2. Neutrosophic observation and terminology 


Numerous kinds of neutrosophic observations were provided, such as quantitative neutrosophic 
data that suggested a number might fall within an unknown interval (a,b), [9]. There are several 


ways to represent the interval value of a neutrosphic number. [16] and [17] have used neutrosophic 
interval values as Z, =Z,+Z,I, , where, I, €[/,,/,,]. In addition, we apply the same notations 


from [17] to the considered neutrosophic data, which are represented in the intervals form as 
Zy €[Z,,Z,], where, 'Z,'is the lower value and ‘Z,,’ is the upper value of the neutrosophic 


variable Z,. We choose a neutrosophic sample of size n, €[n,,n,] using the simple random 
sampling without replacement (SRSWOR) technique from a finite population of 'N' units 
(P,P,,...,.P,). Let yy@eLy,,y,] be the i” sample observation of our neutrosophic data. The 


auxiliary variable x, (i) is similarly represented as x, (i) €[x,,x,]. 

Let y,(@€Ly,,y,] and x,(@ €[x,,x,] be sample mean for the neutrosophic y, and x, 
respectively and neutrosophic auxiliary variable x, is correlated to our neutrosophic study variable 
yy The population means for the neutrosophic variables y, and.x,, which represent the overall 
averages of the neutrosophic data set, are represented by Y, E Y,.¥,] and X ve [xX ox | 
respectively. For y, and x,, the neutrosophic coefficients of variation are C\, €[C,,,,C\yy] and 
Cw €[Cwr»Cwy] respectively. Furthermore, the neutrosophic correlation coefficient between 
neutrosophic variables x, and y, is denoted as /,,y € [ Pry Pyat A The neutrosophic coefficient 


of kurtosis for auxiliary variable x, is denoted by £,..)y €[Py LA. uy]. The neutrosophic mean 


errors are €,, €[€,,,€,,] and &,, €[€,,,€,,]. In order to find the estimator’s bias and MSE, we 


write 

YwO=%y (1+2yv), Fy@O=Xy (14%), 
E(@.y)=E(@,y) =0 and 

E(@yy) = OyCoy > EC) = Py Cov s ECF av) = Oy Pye Cy Con - 


= aj es _ - = 2 2 2 —2 2) 2 
Where, €y €[@,,@y]. Gu €l@n.@u1, ww €1%282>20Sul, Sn Elen eu], Sy ele, eu). 


2 
oO 
2 _ Onn 2 _ ~ yN = 
Cw ~~ xX > Gww ~ y > Cy = PywCywCyn > 


where, Cc. € [C°,, Coy] ’ Ci € Kener Cyl : 
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o N 
= yx 
Pry Px SPP als 


yN Oo 


t= 
0, Unf) where, 0, €[,,9,1, ny €ln,.my], oy €lo,,. oy], Oy €[Oy Oyu]; 
Ny 
Onn € [Ov Ornu | . 
Likewise, in a _ neutrosophic environment, bias and MSE can be defined as 
Bias(y,) €[Bias,, Bias,,], MSE(y,) €[MSE,, MSE,, | respectively. 


The graph below illustrates the flowchart of the proposed study employing neutrosophic numbers. 
The flowchart below has been recommended by [16]. 


Simple Random 
Sampling 


Determinate a 
Data 


3. Some existing neutrosophic estimators 


The estimation of neutrosophic population mean Y, using neutrosophic sample mean which is 


provided by, 

ty =n (3.1) 
The variance of the estimator ¢,,, to the first order of approximation is provided by, 

VG SOY: 9 Whee, ti Eph] (3.2) 
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[16] recommended the usual neutrosophic ratio estimator of neutrosophic population mean Y,, 


utilising the known neutrosophic population mean of X as, 


ty = Yy (2) (3.3) 
The bias term and MSE of the estimator f,,, are given by, respectively 

Bias(t,y) = OyYy (Coy ~Cyav) (3.4) 
MSE( ty) = Oy¥q (Coy + Coy —2C uv ) (3.5) 


Where, Cy.y =PywOwGyy and thy €[ty,.ty]. 


Motivated by [2], the neutrosophic product estimators is given by, 


ty = Vy (=| (3.6) 


The bias and MSE of the estimator's (¢,,,) are provided by, 


Bias(t,y) = Oy¥yC 


yxN (3.7) 
MSE(tyy) = Oy¥q (Coy + Coy +2C yay) 5 Where, thy Ely sty ]- (3.8) 


Inspired by [4], [16] proposed the following neutrosophic ratio estimator that utilises the coefficient 
of variation (CV) of the neutrosophic variable X as, 


At Me ae 
tay =n [8 | (3.9) 
N xN 


The bias and MSE of the estimator t¢,,,, up to the first order of approximation are respectively, 


Bias (47) = Oy¥y AinCoy — Ain Coan) (3.10) 
MSE(t,) = OY. (Cy + re ea a DAAC x) (3.11) 
where, 4, == aS and ty [ty sty]. 


nt xN 


Motivated by [3], [16] suggested the following neutrosophic exponential ratio-type estimator as 
follows: 


Xe 
t., = y, exp| —~— 3.12 
sn — Yn »( | ( ) 


The estimator's (¢;,) bias and MSE are provided by, 
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; = (3 1 
Bias(ts, ) = Oy Yy (3 C.,- a (3.13) 
= C 
MSE(t;,,) = Oy¥x, Ge + ai 2 Cn , where, t,, €[t., ts]. (3.14) 


Motivated by [3], [18] suggested the neutrosophic product-type exponential estimator as, 


ok 
t= exp| ——__* 3.15 
6n — Yn (BF | ( ) 


The bias and MSE of the estimator ¢,,, are given by, 


, aes. 1 
Bias(tey ) = 9y¥y (Fe: + $C] (3.16) 
_ C 
MSE (tgy ) = 9yY¥,, e + a ee > wheres tlt: (3.17) 


Motivated by [19], [17] suggested the neutrosophic ratio estimator ¢,,, utilising the known 


1N? 


population coefficient of correlation p,,, is provided as, 


yxN 


Xo+ 
pn | (3.18) 
Xv + Py 


The bias and MSE of the estimator ¢,, are respectively given as, 


Bias(t,y) = OyY¥y AgyCon — Aw Cyan) (3.19) 
MSE(tyy) = OyYy (Cow + Aan Coy ~ 2a C un) (3.20) 
where, 1,, ==—~* and! @,-eltzgct 5 |. 


Xv + Py 


Inspired by [20], [17] proposed the neutrosophic ratio estimators utilising CV and the coefficient of 
kurtosis of auxiliary variable X, expressed as follows: 


X.+C 
tn =39| Be : 7 (3.21) 


Bacon *n a Cw 


The estimator's (7,,) bias and MSE are defined as, 
Bias(tyy) = Oy¥y (As Coy —AayC puny) (3.22) 


MSE (tay) = Oy¥y (Coy + AswCoy — 2g Cyan ) (3.23) 
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Bacon X w 


where, A,,. = ———=——__ 
BocowX n + Cw 


and ty €[tgrstgy1- 


Motivated by [8], [16] recommended the neutrosophic modified exponential ratio estimator defined 
as, 


= (aX , +b) — (ax, +b) 
= Ree Ne A 24 
ead (aX , +b) +(aXy +B) om 


where a and b denotes the auxiliary neutrosophic parameters. 


The bias and MSE of the neutrosophic exponential ratio estimator f,,, are respectively given as, 


23 1 
Bias(tyy ) = OyYy F hayCoy ~ 5 seus (3.25) 
MSE(isy) = Oy¥y (Coy + Ay Coy — 2a Cen } (3.26) 
aX 
where, A,, * Fak, +b arid 1, 5:€ [tgp Stay 


From (3.22) and (3.25), we see that the estimators ¢,, and f,, recommended by [17], [16] 


respectively, are biased estimator. Bias is disadvantageous in certain applications. 


Motivated by [8] and [17], we have suggested an almost unbiased estimator for the population mean 
Y,, using neutrosophic auxiliary information in SRSWOR. 


4. The proposed neutrosophic almost unbiased estimator 


— K aS - 
Suppose t) =Vy 5 t=] Vy | , = , oo[ (SF +b)—(axy 2) 


ax, +b (aX , +b) + (ax, +b) 


such that f),t,,t,¢A, where Ais the collection of all possible estimators for estimating the 


neutrosophic population mean Y,,. Justification of using a and b in proposing estimators is given 


by [22]. By definition, the set A is a linear variety if 


2; 
t=) at,<A. (4.1) 
i=0 


2 
for ya, =l1 anda,eR. (4.2) 


i=0 


where @,(i=0,1,2) refers to the statistical constants, R stands for the set of real numbers and K 


takes distinct values, that is, K = 1 for the ratio estimators and K =~-—1 for the product estimators. 
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The bias and MSE of the estimators ¢, and t, are, respectively, given by (up to the first order of 


approximation) 
~| 5-9 { K(K+1) 
B(t,) = Oy¥y | A Cow an —AK Py Cu (4.3) 
MSE(t,) = Oy¥q | Coy + A?K Coy — 24K PyuyCyyCuy | (4.4) 
= 2K+K*\ AK 
B(t,) =O,Y,, [cs P| — AE pELu (4.5) 
= A’ K* 
MSE(t,) = 0,Yy [cx - Cy - RPC Cu (4.6) 
Ax 2 
where, 2 =—= ; a and b denotes the auxiliary neutrosophic parameters. 
aX, +b 


The bias and MSE of the estimator f are determined using, 
Vy) =%y (1+2yO), Xy@=Xy (1+Fy@). 
Such that, E(é@,y)=E(@,,)=0. 


Defining the estimator ft in terms of 2, and é,y » we obtain, 
X, +b) 
. = ze = = aX y+ 
=a,Yn(1 vam = ay 
faao( v5} (ty) | | 
2 X,+b)—(aX,(1+2@,)+b 
+, Vy (1+2,y exp K caste Coa Ce) (aX Yeon) 
: (aX , +b) +(aXy (1+ )+b) 
= a,¥v (1+2,y) +a] Pw (1+2y)(1+ 2%)" | 


7 “1 (4.7) 
+a,Ywn (1+ ex}os| K{ 45, [1+ 42. ] 


ie 
aX, +b 


where, 1 = 


Expanding the terms on the right side of the equation (4.7) and having the terms up to the first order 
of approximation. 
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(a, +a, +a, )+(a, +a, +a, Je -1Ké., [a + | tale,” 
f=y5 , (4.8) 
aoa aus aAe.x ew i 2 ANE Ey 


K(K a 
2 


The bias of the estimator ft’, up to the second power of @'s , can be computed by subtracting Yv 
from both sides of equation (4.8) after taking the expectation of each side. 


* i : K a 
B(t’) = 0, ¥vC2, aa aa +a,0 PR) -{ ka,A + Sash Pros (4.9) 


xN 


From equation (4.8), we have 
-Yv =¥v| Oy —AaKe,y | (4.10) 
a, 
where, a@ =| a, ser ‘ (4.11) 
By taking square of each side of (4.10) and then taking expectations, we get MSE of the estimator 
t’, to the first order of approximation as, 
MSE(t')=0,Yw[ C2, +C2,AKa(AKa-26, ) | (4.12) 


Cw 
where, dy = Pyxw a 


xN 


Minimum MSE of the estimator t would occur when 


pan (4.13) 


N 


Applying this value of @ ne in equation (4.1), we get the optimum value of estimator 1 
(optimum). 
Hence, the minimum MSE of the estimator f is provided by 


min .MSE(t’)=0,¥wC?, (1-p2y) (4.14) 


This MSE is equivalent to the classical linear regression estimator. 


By using equation (4.11) and (4.13) we have 


by ase Os (4.15) 
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There are only two equations with three unknowns in equations (4.2) and (4.15). The unique 
values for as, ({=0,1,2) can't be found. In order to get unique values for as , we will apply the 


linear constraint 


5 a,B(t,) =0. (4.16) 


sth 


where B(t,) represents the bias of i" estimator. 


The matrix form of equations (4.2), (4.15), and (4.16) is as follows: 


Ge 1 
0 1 5 jla|qle (4.17) 
0 Bt) Bit, |) \° 
where, a@=—_. 

AK 


By using equation (4.17), The unique values of as ((=0,1,2) are obtained as, 


B(t,)- ne )_ @B(t,)+aB(t) 
(04 = 
a _ Bu) 
ra 
aB(t,) 
2 (4.18) 
Bit) 
Be) 
aoe —aBtt,) 
| Bu) 
BU) ~~ | 
such that, 


2 
Yia,=1 and a,eR. 


i=0 
By using these a's, (i = 0,1, 2) eliminate the bias up to terms of order o(n'') at (4.1). 


5, Some member of the suggested estimator 't° 


Several estimators of the population mean are provided by choosing the appropriate combination of 


constants a, b and as ,(i =0,1,2), as shown in the following scheme. 


Estimator Values of 
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ay 1 2 a b K 
try =v 1 0 0 0 0 1 
The mean per unit estimator 
mee 0 d 0 1 0 1 
Ly =Yn| a 
Xy 
[16] estimator 
mee 0 1 0 1 0 i 
bby = Yu| S 
X y 
Neutrosophic product estimator 
_ {Xy+C 0 1 0 1 C, 1 
lay = Yu = = : 
coe a Oe 
[16] estimator 
X._-< 0 0 1 1 0 1 
te = Vy exp| —“——* 
sn — Yn »(2 =a 
ee 0 0 1 r 0 4 
toy = Yy exp| ~——* 
6N — YN (= ee) 
[18] estimator 
we ie Pie 0 1 0 1 poet 
7N N = 
Xy + Py 
t a) Bese af Cw 0 1 0 Pian Cw 1 
8N = 
i =i Pantn a Cw 
[17] estimator 
X +b)—(ax,, +b 0 0 1 a b 1 
toy = Vy EXP iC =" ee) 
(aX , +b) +(axy +b) 


In addition, the estimators mentioned above, by simply changing the values of a, b and as( 


i=0,1,2). From the proposed estimator f° at (4.1), several estimators can be created. 


6. Empirical Study 


Population 1 : Here, we use the data provided by [21], which is birth rate and data on the natural 
growth rate obtained from the sample registration system (SRS)(2020) is used to demonstrate the 
performance of different estimators of neutrosophic study variable Y, €[Y,.Y,,]. The neutrosophic 


auxiliary variable and neutrosophic study variable are given by, 


X, : Birth rate Y, : Natural growth rate 
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The PRE formula is given by 


MSE(estimator) 

Table 1. Descriptive statistics of birth rate and natural growth rate data as per SRS 2020. 
Parameter Neutrosophic Value | Parameter Neutrosophic Value 
Ny [36, 36] Nig [10, 10] 

Low [14.472, 18.053] Co [0.333, 0.411] 
fx [10.894, 10.631] ne [0.620, 0.849] 
Cx [0.213, 0.236] ae [3.034, 2.173] 
Table 2. Neutrosophic MSEs and PRE of the different estimators. 
Estimator MSE I PRE 
ae [0.951886, 1.378714] | [0, 0.31] [100, 100] 
bs [0.490003, 0.586746] | [0, 0.16] [194.26, 234.98] 
a [2.092575, 3.174567] | [0, 0.34] [43.43, 45.49] 
bis [2.070519, 3.145647] | [0, 0.34] [43.83, 45.97] 
tei [0.672372, 0.820965] | [0, 0.18] [141.57, 167.94] 
i [1.425287, 2.163247] | [0, 0.34] [63.73, 66.79] 
te [0.510456, 0.586471] | [0, 0.13] [186.48, 235.09] 
bis [0.492618, 0.586646] | [0, 0.16] [193.23, 235.02] 
ty (a=1,b=1) 0.684577, 0.844600] | [0, 0.19] [139.05, 163.24] 
t [0.385639, 0.586416] | [0, 0.34] [235.11, 246.83] 


Population 2 : Here, we use the data sets considered by [18]. Numerous data points, such as the 
minimum, maximum and mean temperatures over a wide time duration, are included in this dataset 
regarding the All India Seasonal and Annual Temperature Series. The neutrosophic auxiliary 
variable and neutrosophic study variable are given by, 


X, : Minimum and maximum temperatures in January and February. 
Y, : Temperatures from March through May. 


Table 3. Descriptive statistics of the population 2. 


Parameter Neutrosophic Value | Parameter Neutrosophic Value 
Ny [119, 119] he [20, 20] 

im [13.908,24.674] Cy [0.0007, 00.0006] 
x [20.694, 31.560] Pao [0.6274, 0.7202] 
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[0.0017, 0.0015] Bacon [4.973, 5.583] 


xN 


Table 4. Neutrosophic MSEs and PRE of the different estimators. 


Estimator MSE ti PRE 

fs [8.87822E-06, 1.72193E-05] | [0, 0.48] [100.00, 100.00] 
Be [3.39108E-05, 5.58304E-05] | [0, 0.39] [26.18, 30.84] 

on [8.78263E-05, 0.000173977] | [0, 0.50] [10.11, 9.90] 

a [3.39013E-05, 5.58219E-05] | [0, 0.39] [26.19, 30.85] 

i [8.39692E-06, 1.21037E-05] | [0, 0.31] [105.73, 142.26] 
es [3.53547E-05, 7.1177E-05] | [0, 0.50] [25.11, 24.19] 
ie [3.06832E-05, 5.20435E-05] | [0, 0.41] [28.94, 33.09] 
i [3.39089E-05, 5.58288E-05] | [0, 0.39] [26.18, 30.84] 
ty (€=1,b=1) | [7.61584E-06, 1.13888E-05] | [0, 0.33] [116.58, 151.19] 
a [5.38372E-06, 8.28831E-06] | [9, 9.35] [164.91, 207.75] 


7. Discussion 


In Table 2, as can be shown that the sampling variance of the neutrosophic estimator ¢,, 1s 
[0.951886, 1.378714] and the neutrosophic exponential ratio estimator’s (1,,) MSE is [0.672372, 
0.820965] whereas the range of the neutrosophic MSEs (t,,,) is [2.070519, 3.145647]. The reason 
behind the higher MSEs of the neutrosophic estimators t,,, over the neutrosophic estimator f,, is 


that the low value of the neutrosophic CV of the neutrosophic auxiliary variable x. In comparison 
to the other population mean (Y,,) estimators in competition, the MSE of the proposed class of 


neutrosophic estimators (t’) is the lowest at [0.385639, 0.586416] with higher PRE. 


Similarly, the proposed and existing estimators' MSEs and PRE are given by Table 4. Additionally, 
it has shown that the suggested estimators' MSEs and PREs are lower and higher respectively, in 
comparison to those of the other estimators already in use. Its value is highlighted by the bolded 
text. Hence, Tables 2 and 4 confirmed these results, showing that the suggested estimators are 
preferable to the existing ones. 


8. Conclusion 


In this paper, we have recommended an almost unbiased neutrosophic estimator to estimate the 
neutrosophic population mean Y, y using some known value of auxiliary information. In order to 


evaluate their precision, we computed bias and MSE for the suggested estimators, concentrating on 
first-order approximations. We evaluated our suggested estimator with existing estimators using 
data from a natural environment. Via conducting empirical study, we have identified strong 
evidence that our suggested estimator performs better than the existing ones within the context of 
neutrosophic sampling. It should be emphasised that in cases where the study variable's observations 
are nondeterministic or uncertain, the neutrosophic estimators are best suited for enhanced 
population mean estimate. Therefore, it might be advised to use the proposed class of estimators to 


Rajesh Singh, Anamika Kumari, Florentin Smarandache and Shobh Nath Tiwari, Construction Of Almost Unbiased 
Estimator For Population Mean Using Neutrosophic Information 


Neutrosophic Sets and Systems, Vol. 76, 2025 462 


estimate the neutrosophic population mean Y,, in a variety of real-world settings, including those in 


the social sciences, economics, agriculture, biology and mathematics. 
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